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r^j In this work we compute the color fields in the mediator plane between a 

static quark and a static antiquark using quenched lattice QCD. In special, 
we see the effect of the quark-antiquark distance on the flux tube. To 
CO obtain this results an improved multihit technique is developed and an 

extend smearing technique is used. Then, we also discuss the flux-tubes 
[ in a system composed of two quarks and two antiquarks. The ground and 

+^ first excited states fields are studied for different dispositions of the system. 

PACS numbers: 12.38.Gc, 11.15.Ha 

r^, 1. Introduction 

I This work is divided in two parts. In the first one we will study an 

^ isolated fundamental flux-tube in a quark-antiquark system. In the second 

one we will study the fields a system made of two quarks and two antiquarks, 
where we can observe the interaction between two fundamental flux tubes. 

qq 2. Chromo-fields computation 

^ To calculate the gauge invariant squared chromoelectric and chromomag- 

netic fields in the lattice, we only have to calculate the correlation of the 
Wilson loop operator W which in large euclidean time limit t —> oo should 
describe the static system, with the plaquettes corresponding to each one 
of the fields. Concretely, we have: 



(E?> = <%> " ^ (1) 

(Bl) = <VLM- (Pil) (2) 
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where = 1 - ^Tr[l7 /x (s)J7 I/ (s + v)ul(s + u)ut{s)]. The spatial indices 
j and /c complement index i. The Lagrangian density is given by C — 
\{E 2 -B 2 ). 

3. Fundamental flux-tube 

For a simple, quark-antiquark singlet, the Wilson Loop operator W de- 
scribes the system. We compute the chromo-fields in the mediator plane 
between the quark and the antiquark, using 1100 pure gauge 32 4 configu- 
rations with /3 = 6.0 . In order to improve the signal to noise ratio, we use 
two techniques: an improved version of the multihit [6] and an extended 
spatial smearing technique. 

3.1. Signal to noise ratio improvement 
In the multihit method we replace each temporal link by its thermal 
average, with it's first neighbors fixed, that is Ua Ua = - dU4 U4 % — x — . 

5 ' 5 ,44 f dU4 e ?Tr [U 4 Ft] 

We generalize this method by replacing each temporal link by it's ther- 
mal average with the N th neighbors fixed, that is: 

In this way we have an error reduction greater than the one of multihit. 
To increase the ground state overlap, we use a spatial extended smearing 

Ui v su(3) 5j- +^ 2 ^ s ij 4] ( 4 ) 

3 j 3 

with staples S]p Sf- and Sfj being given on figure [TJ 

Even using this technique we weren't able to find a value of t for which 
the plaquette to Wilson Loop correlators are stable within error bars, while 
still have a sufficiently high signal to noise ratio. To solve this, we note that 
the correlator which gives the average of field (F) should be given by the 
formula (F) t = (F)^ + be~ At for large values of t, with A = V\ — Vb, being 
the difference between the first excited and the ground state potentials. To 
compute A, we a use a variational basis of four levels of APE smearing, 
with the potentials V\ and Vq being given by the solution of the variational 
generalized eigensystem 

(Wij{t))<2{t) = ™ n (0<Wy(0)> C ?(t) (5) 
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Figure 1. Staples used in the extended spatial smearing. 




Figure 2. Left: Lagrangian density in the mediator plane between the quark and 
the antiquark. Right: Width of the flux tube W\/ e as a function of R. 



where (Wij) = (Oi(t)O^(0)) is the correlation between the meson anni- 
hilation and creation operators at time t and in the smeared states i and 
j respectively. 



3.2. Results and conclusions 

The results for the Lagrangian density in the quark- antiquark mediator 
plane are shown on the left side of Fig. [2j As can be seen the tube flux 
becomes wider as the quark-antiquark distance is increased. To estimate 
the width of the flux tube we define the value Wi/ e by (C)(wi/ e ) — ^ (C)(0). 
The results for this are shown on the right side of Fig. [2| 

It was predicted [U ?] that the squared width w 2 of this flux tube di- 
verges logarithmicaly as R 00, that is w 2 ~ Wq log , where w 2 . This 
behavior is called "roughening". So, we can say that the flux tube rough- 
ening has been observed from R ~ 0.4 fm to R ~ 1.2 fm. However, we 
weren't yet able to confirm the that the behavior is logarithmic. 
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Figure 3. Left: Anti-parallel geometry. Right: Parallel geometry. 



4. QQQQ system 

Now, we consider a system of two quarks Q\ Q2 and two antiquarks 
Q3Q4. The study of this kind of systems is important for the under- 
standing of meson meson scattering processes and the possible formation 
of tetraquarks, particles made of two valence quarks and two valence anti- 
quarks. 

In this system we have two linearly independent color singlets. Those 
can be, for instance, the two meson states: \I) — ^IQiQjQ^j) and |77) = 
\\QiQjQjQi), or the anti-symmetric and symmetric color states: \A) — 

— \H)) and \S) = + \II)). In accordance with lattice results 

PQ[5], the ground state could be either 1 77) or |A), with the static poten- 
tial being given by the flip-flop potential Vff = min(V/, Vn, Vr), where Vi 
and Vn are the two possible two-meson potentials, given by the sum of the 
intra-meson potentials, and Vt is the tetraquark potential, which confines 
all the particles, with the confining part being proportional to the minimal 
length of a fundamental string linking the four particles. 



Geometries 

We use the two geometries, shown in Fig. [SJ In both, the four particles 
form a rectangle, however in one of them — the parallel — similar particles 
are on the same side of the rectangle while on the other they are on opposite 
corners — the anti-parallel. With the first one, we can study the tetraquark 
to mesons transition while with the second we can see the transition between 
the two meson states. 

In this work, we use 1121 Quenched 24 3 x 48 lattice configurations with 
f3 = 6.2. APE Smearing was applied to the spatial links, while Hyper-cubic 
blocking [3J was applied to the temporal links. 
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Figure 4. Left: Loops used to study the antiparallel geometry. Right: Loops used 
to study the parallel geometry. 




Lforr t = 6 and r 2 = 8 Lforr t = 6 and r 2 = 10 




Figure 5. Lagrangian density for the ground state of the antiparallel geometry. 



4-2. Variational method 

To obtain not only the ground state, but also the first excited state, we 
use a variational basis similar to Eq. ^ but using two different kind of 
operators. This way, the Wilson loop which appears on Eqs. and ^ 
is W n = c l n WijCn. The two operator basis will be formed by the \I) and 
\II) annihilation operators in the anti-parallel geometry and in the parallel 
geometry, by the \I) and \A) [U [5] annihilation operators. This gives the 
matrix elements on Fig. [4| 



4.3. Results 

The results for the Lagrangian density C in the ground state and the first 
excited state are given on Figs. [5] and [6| For the ground state the system 
collapses in a two meson state, when T2 ^> n., as expected. Looking at the 
Wilson loop composition we also conclude that when r\ — T2 the ground 
state is color symmetric. The excited state is not so readily explainable. 
The results for C in the parallel geometry can be seen on Figs. [7] and |8j 
The results are the expected ones for the ground state, with the system 
passing from a two meson state to a tetraquark state as we increase V2> We 
discuss the first excited state in the conclusion. 



4-4- Discussion/Conclusion 

The results for the ground state are consistent with the flip-flop ground 
state potential with the color fields dispositions being consistent with 
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Lforr, = 6andr 2 = 6 L for r- = 6 and r 2 = 8 Lforr, = 6 and r 2 = 10 




Figure 6. Lagrangian density for the first excited state of the antiparallel geometry. 



Lforr, =8 and r 2 = 8 L for ^ = 8 and r 2 = 10 Lfor r t = Sand r 2 = 12 Lforr, =8 and r 2 = 14 




Figure 7. Lagrangian density for the ground state of the parallel geometry. 



\II) or \A) where expected. 

We can explain the excited states if we consider that the vibrational and 
rotational flux tube excitations effects are negligible. In this case the first 
excited state is orthogonal to the ground state. In that case, if we compute 
the Casimir factors for this orthogonal states, as in [2], that they predict 
the correct behavior of the fields, with a repulsion between particles, where 
we expect it. 
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Figure 8. Lagrangian density for the first excited state of the parallel geometry. 
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